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I. ABSTRACT 


It is important that accurate estimates of crew exposure to radiation are obtained for 
future long-term space missions. Presently, several space radiation transport codes ex- 
ist to predict the radiation environment, all of which take as input particle interaction 
cross sections that describe the nuclear interactions between the particles and the shield- 
ing material. The space radiation transport code HZETRN uses the nuclear fragmentation 
model NUCFRG2 to calculate Electromagnetic Dissociation (EMD) cross sections. Cur- 
rently, NUCFRG2 employs energy independent branching ratios to calculate these cross 
sections. Using Weisskopf-Ewing (WE) theory to calculate branching ratios, however, is 
more advantageous than the method currently employed in NUCFRG2. The WE theory 
can calculate not only neutron and proton emission, as in the energy independent branching 
ratio formalism used in NUCFRG2, but also deuteron, triton, helion, and alpha particle 
emission. These particles can contribute significantly to total exposure estimates. In this 
work, photonuclear cross sections are calculated using WE theory and the energy indepen- 
dent branching ratios used in NUCFRG2 and then compared to experimental data. It is 
found that the WE theory gives comparable, but mainly better agreement with data than 
the energy independent branching ratio. Furthermore, EMD cross sections for single neu- 
tron, proton, and alpha particle removal are calculated using WE theory and an energy 
independent branching ratio used in NUCFRG2 and compared to experimental data. 

II. INTRODUCTION 

Radiation protection is an important technology for the future of space exploration. When 
traveling into space, radiation is emitted from a variety of sources, such as the Van Allen 
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trapped radiation belts, solar energetic particles, or galactic cosmic rays. Accurate estimates 
of the radiation exposure to astronauts or electronics is important for long term missions, 
such as the exploration of Mars. Several particle transport codes have been developed that 
estimate space radiation exposure. These transport codes take as input particle interaction 
cross sections that describe the nuclear interactions between the particles and the shield- 
ing material. Therefore, an accurate and efficient method of calculating cross sections is 
fundamental for shielding optimization. 

A relativistic nucleus-nucleus collision will result when a cosmic ray nucleus interacts with 
a spacecraft wall or shielding material. Nucleus- nucleus collisions can be mediated by either 
the strong or electromagnetic (EM) forces. EM processes dominate when the collision’s 
impact parameter is larger than the range of the nuclear force, so no nuclear interaction 
occurs. A reaction proceeding via the EM force is often called Electromagnetic Dissociation 
(EMD). For this type of reaction, a projectile nucleus is excited by the absorption of a 
virtual photon, which is supplied by the target nucleus. The projectile subsequently decays 
by the emission of one or more nucleons, with neutron emission the primary decay mode for 
nuclei with high atomic numbers, Z, and for photon energies below a few tens of MeV. In 
the region of low Z and energies above a few tens of MeV, light ion production will compete 
with neutron emission. The Coulomb barrier suppresses light charged particle emission for 
heavy nuclei, but for lighter nuclei, the Coulomb barrier is lower, thus allowing for light ion 
production. 

The high charge and energy (HZE) particle transport code HZETRN is a computational 
tool used for space radiation studies and shield design. HZETRN uses a deterministic ap- 
proach for particle transport to enable fast and accurate estimates of the relevant dosimetric 
quantities [1, 2], HZETRN uses the nuclear fragmentation model NUCFRG2 [3] to calcu- 
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late EMD cross sections. Currently, NUCFRG2 employs an energy independent branching 
ratio for the calculation of EMD nucleus-nucleus cross sections. This energy independent 
branching ratio can only account for neutron and proton emission. In this paper, a way 
of calculating EMD nucleus-nucleus cross sections, using an energy dependent branching 
ratio, is presented. This branching ratio can account for the emission of a neutron, proton, 
deuteron, triton, helion, and alpha particle, which is a significant improvement over the 
current NUCRFG2 branching ratio model. It should be noted that the additional ions have 
been shown to account for a non-negligible portion of the total exposure in some cases. The 
energy dependent branching ratio, which is calculated from Weisskopf- Ewing (WE) theory, 
is first folded into the photonuclear cross section formula and then incorporated into the 
EMD nucleus-nucleus cross section equation. 

A large part of the uncertainty in an EMD nucleus-nucleus cross section calculation comes 
from the photonuclear cross section. The uncertainties generally arise due to the poor the- 
oretical treatment of the branching ratio. This paper will compare an energy dependent 
branching ratio, which is calculated from WE theory, to the energy independent branch- 
ing ratio currently employed by NUCFRG2. It is shown here that the energy dependent 
branching ratio, calculated with WE theory, provides better estimates of the photonuclear 
cross section, which consequently give more accurate EMD nucleus-nucleus cross sections. 
A review of the photonuclear cross section is first given, followed by a summary of the 
WE method and the energy independent branching ratio used in NUCFRG2. The EMD 
nucleus-nucleus cross section is also reviewed and comparisons are made. 
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III. PHOTONUCLEAR CROSS SECTION 


This paper focuses on photonuclear reactions where a photon interacts with a nucleus, 
typically at an energy above a few MeV. When a nuclear excitation occurs, the main contri- 
bution comes from the Giant Dipole Resonance (GDR); the reaction results in the emission 
of neutrons, charged particles, or gamma rays. The photonuclear cross section is a funda- 
mental component of an EMD nucleus-nucleus cross section, as it describes the response of 
the excited nucleus to spectator photons. If there is uncertainty in the photonuclear cross 
section, it will ultimately induce uncertainty into the EMD cross section. Therefore, it is 
important that the photonuclear cross section is well understood. 

The photoparticlc total cross section for producing particle b is given by the product of 
the photonuclear absorption cross section, cr^bsiE^), and the branching ratio, gj , , representing 
the probability of a specific decay mode [4, 5]; 


*^"(■^7 ) E) gb ^"abs(-Ey) > 


( 1 ) 


where E 7 is the photon energy. 

The branching ratio, gb, will be discussed later, and it will be shown that the branching 
ratios can be modeled as either energy independent or energy dependent. Note that the 
branching ratio in equation (1) is energy independent; the dependence will be explicitly 
written if it is included. 

Typically, experimental data are used for the absorption cross section, when available. 
For reactions with no experimental absorption cross section data, parameterizations are 
often used [5]. In addition, applications which require rapid cross section calculations, such 
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as space radiation transport codes, employ a parameterized absorption cross section, since it 
would be impractical to put experimental photonuclear cross sections into EMD calculations 
for transport code analyses. The parameterized absorption cross section used in this paper 
will now be discussed. 

A. Absorption Cross Section 

The photonuclear absorption cross section, cr a b s (-E 7 ),is parameterized in the region near 
the giant dipole resonance as 


a ‘ b ’ ( ' El) 1 + [(£? - ££ DR )7Ey 2 ] ' (2) 

The abbreviation, GDR, stands for the giant dipole resonance, as mentioned previously. 
Here, E GDR is the energy at which the photonuclear cross section has its peak value, and 
T is the width of the electric dipole (El) giant dipole resonance. Values of T used in this 
paper can be found in Table I. Also, 
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with the Thomas-Reiche-Kuhn cross section given by [5] 
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with the * referring to the excited nucleus. Here, A /"*, Z*, and A* are the neutron number, 
atomic number, and the mass number of the excited nucleus, respectively. The GDR energy 
is [5] 


with 
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and 


R 0 = r o Al * 1/3 . 


(7) 


The parameters are [5]: 
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e = 0.0768 , (8) 

Q' = 17 MeV , (9) 

J = 36.8 MeV , (10) 

ro = 1.18 fm , (11) 

177 ' 0.7 = 0-7 lUnucleon 5 ( 12 ) 


where m nuc i eon is the nucleon mass, taken as 938.95 MeV. 

IV. BRANCHING RATIO 

A branching ratio provides the relative probability of a specific type of decay. In this 
section, a form for the branching ratio, which is dependent on the excitation energy of the 
compound nucleus, will be introduced that is calculated on the basis of the WE evaporation 
model of nuclear reactions. In addition, the energy independent branching ratio, used in 
NUCFRG2, will be presented. This energy independent branching ratio is not restricted to 
compound nucleus decay, but can also be applied to direct reactions. 

A. Energy Dependent Branching Ratio 

Weisskopf- Ewing (WE) theory is a commonly used method for computing an energy 
dependent branching ratio of a compound nucleus reaction. In a compound nucleus reaction, 
a compound nucleus is formed in an excited energy state and then decays into the products of 
the reaction once reaching statistical equilibrium. The formation and decay of the compound 
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nucleus can be considered independent processes because only energy, angular momentum, 
and parity determine the de-excitation of the compound system. 

The Hauser-Feshbach (HF) method provides another way of calculating cross sections 
for compound nucleus reactions. The difference between WE and HF lies in the number 
of allowable dynamical paths. The HF method accounts for more states, since parity and 
angular momentum are considered. The WE theory provides a simple way of calculating 
compound nuclear reactions going to continuum states, because only energy is conserved. 
The HF method has greater complexity than the WE theory and presents a much greater 
computational cost. Despite the fact that the HF method is generally favored for calculations 
of reactions proceeding through a compound nucleus state, the WE theory is chosen here 
due to its straightforward formulation and relatively low computational cost. Thus, the 
WE theory can be seen as a good alternative to the HF method for applications, such as 
radiation transport codes, that require near real time simulations. Here, computational 
efficiency must be weighed against the sensitivity of the results to the model. 

The WE theory is developed by first considering a nuclear reaction where a particle a 
strikes a nucleus A to produce a residual nucleus B and an outgoing particle b. The reaction 
proceeds though a well-defined compound state, 


a + A -»• C* -»• B + b , 
(State A) (State C ) (State B) 


( 13 ) 


where C represents the intermediate compound nucleus and the * signifies an excited state. 
Particles a and b may be particles, such as photons, neutrons, and protons, although they 



could also be nuclei, such as deuterons, tritons, helions, or alpha particles. Note that a 
photonuclear reaction proceeds via r j + A^ A* ^B + b. 

The probability that the compound state will decay by emission of particle b is given by 
the energy dependent branching ratio [6] as 


gt(Ei) 


rygc) 
r«(£y ’ 


(14) 


where T& is the partial decay width for emitting particle b, T to t is the total width of the 
compound state, and E £ is the excitation energy of the compound nucleus C. It is important 
to note that for a photonuclear reaction, the excitation energy will correspond to the photon 
energy, E~. f . Equation (14) will be referred to as the WE branching ratio. The total width is 
defined as the sum of all the partial widths. This work considers the emission of six particle 
types: neutron, proton, deuteron, triton, helion, and alpha. Emissions of other particles are 
assumed negligible. As a result, the total width is given by the sum of the six individual 
decay widths 


r«(sy = r „(££) + T r (E- c ) + r„(sy + r,(Bj) + r h (E’ c ) + r „(sy , (is) 

where n, p, d, t, h, and a represent a neutron, proton, deuteron, triton, helion, and alpha 
particle, respectively. Other decay modes, like fission and gamma emission, could also 
be included in equation (15) as additive terms. These decay modes have been discussed 
elsewhere [7, 8] and will not be discussed here. 
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According to WE theory, the partial width of the decay in the numerator of equation 
(14) will be contingent on whether particle b is charged. For the neutron, the partial width 
of decay is expressed as [9, 10] 
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(16) 


where M n , s n , and e n are the mass, spin, and binding energy for the emitted neutron, 
respectively. Here, (3 n = (2.12A B 2 ^ 3 — 0.050)/C n and C n = 0.76 + 2.2 , where Ab is the 
atomic mass of the residual nucleus [11]. The geometrical cross section cr geo , is represented 
by the classical target area, cr ge o = 7r(r 0 A B 3 ) 2 , where the radius parameter r 0 is taken to be 
1.18 fin [11], Also contained in equation (16) are the nuclear densities pc and ps, for the 
compound nucleus C and the residual nucleus B, respectively. See references [12, 13] for 
detailed information of all these terms. 

The partial width for a charged particle is defined by [9, 10] 
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where M c , s c , V c , and e c , are the mass, spin, Coulomb barrier and binding energy for one of 
the charged particles ( p , d, f, h, a), respectively. 

Both of the expressions for the partial width contain Tb, the temperature of the residual 
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nucleus after b emission, which is calculated by [9, 10] 


T b = 


( dlnp B \ 1 

V dE* c J 

f E* c -e b \ 1/2 / SEq }B a B 

V Ub(Eq — 6b) / \ E^ — €b a B (Ec — €b) 

X [exp (-7 b[E*c - e b }) x (1 + 7 B [E* C - e b ]) - 1] 


-1 


(18) 


Here, a B , a Bl 5E 0 ^ B , and 7 B are the level density parameter, asymptotic level density param- 
eter, shell correction energy, and damping parameter of the residual nucleus B, respectively. 
These terms and a complete derivation of the neutron and charged particle partial width 
formula, equations (16) and (17), have been discussed in detail in references [12, 13]. 


B. Energy Independent Branching Ratio 

In the previous section, a method was introduced for calculating branching ratios, based 
on WE theory. This branching ratio formalism differs from the method used in NUCFRG2. 
Currently, NUCFRG2 employs a branching ratio from references [14, 15]. This branching 
ratio neglects energy dependence, unlike the WE theory, which takes into account the exci- 
tation energy (or for a photonuclear reaction, the photon energy). This energy independent 
branching ratio is not restricted to compound nucleus decay, but can also be applied to direct 
reactions. A direct reaction has a short interaction time, unlike the long lived intermediary 
state of a compound nucleus reaction. A compound nucleus reaction follows equation (13), 
but a direct reaction proceeds via 
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a ~b A 


B + b. 


(19) 


(State A) (State B ) 

The energy dependent proton branching ratio utilized in NUCFRG2 is given by [3, 14, 15] 


0.5 
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6<Z* A < 8 
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8<Z* a <14 

min 1-95 exp(— 0.075 Z* A ) 

14 < 


( 20 ) 


Here, Z* A and A* a are the atomic number and mass number of the excited nucleus A (either 
the projectile or target), respectively. Due to the dependence on only the atomic number 
and mass number of the excited nucleus A, this energy independent branching ratio can be 
applied to a direct or compound nucleus reaction. 

With the ratio of proton emission now formulated, the neutron branching ratio can be 
found easily from 


9n 1 9p > ( 21 ) 

since only proton and neutron decay are considered in NUCFRG2 [3]. For this crude ap- 
proximation, emission of other particles was thought to compete insignificantly with proton 
and neutron emission, so other channels were neglected for simplicity. It should be noted 
that the above relations assume only single nucleon emission. 


12 



V. NUCLEUS-NUCLEUS CROSS SECTION 


In a nucleus-nucleus collision mediated by the EM force, the target (or projectile) rep- 
resents a source of virtual photons, which impinge upon the projectile (or target). The 
spectrum of virtual photons contains a variety of energies, in contrast to a photonuclear 
reaction where the incoming photon possesses only a single energy. 

The total cross section for electromagnetic nucleus-nucleus reactions can be written in 
the form 


vaa = I dE 7 N(E 1 ) ct 7j4 (.E 7 ) , (22) 

where N(E y ) is the Weizsacker-Williams virtual photon spectrum and a 7y i(E 7 ) is the pho- 
tonuclear total cross section. The Weizsacker-Williams virtual photon spectrum will be 
discussed in the subsequent section. 

Equation (1) expresses a 7 ^(E 7 ) as the branching ratio, g x , multiplied by the absorption 
cross section, (J^JE^,). When replacing the photonuclear total cross section in equation (22) 
with equation (1), the form of the total cross section for an electromagnetic nucleus-nucleus 
reaction will depend on whether the branching ratio is energy dependent or independent. 
For the case of an energy dependent branching ratio, 


&AA 


dE 1 N(E 1 ) gx(E 7 ) cr, b ,(E 7 ) , 


( 23 ) 


whereas an energy independent branching ratio can be brought outside the integral, 


13 



&AA = gx 


(24) 


dE-f N(E y ) o- aba (£^ 7 ) , 


which allows for faster computation. 

A. Weizsacker-Williams Virtual Photon Spectrum 

Equations (22), (23), and (24) require the Weizsacker-Williams virtual photon spectrum. 
The virtual photon spectrum describes the number of photons the target contributes at a 
specific energy and is given by [16, 17] 


N(E 7 ) 


1 2 2 1 
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where Z$ is the nuclear charge of the spectator nucleus (either the projectile or target) and 
Q(f sc is the EM fine structure constant given by CKf sc = e 2 /he. The relativistic beta factor is 
expressed as 



and the relativistic gamma factor of the excited nucleus is given by 


(26) 
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(27) 


1 + 2"lab 

^nucleon 


where 7j a b is the kinetic energy per particle N of the projectile. Equation (25) also con- 
tains the modified Bessel functions, Ko(0 and AT(£), which are functions of the adiabacity 
parameter £; 


-Ay ^rilill 

7' P (He) ' 


(28) 


The minimum impact parameter, £ min , is the value below which the strong force is assumed 
to dominate [18] and is found from the relation [5] 


hriin Ro.l,P T Rq.1,T d over ] a p , (29) 

where Rq.\,p and Rq.i,t are the 10% charge density radii for the projectile and target nuclei, 
respectively. Values of the 10% charge density radii used in this paper can be found in Table 
II. The overlap distance is expressed as ^overlap arid is treated as an arbitrary parameter. 

VI. COMPARISON TO EXPERIMENT 

The primary focus of this work was to evaluate and compare the WE branching ratio to 
the energy independent branching ratio, the method currently used in NUCFRG2, for the 
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purposes of calculating EMD nucleus-nucleus cross sections. To obtain a true comparison of 
these branching ratios, the photonnclear cross section must be examined for reasons which 
will now be explained. Following this explanation, comparisons to experiment will be made 
for both branching ratio methods for EMD nucleus-nucleus cross sections. 

A. Photonuclear Cross Sections using WE and energy independent branching 
ratios 

Photonuclear cross sections, calculated with equation (1) are presented in this section. 
Comparisons to experimental data are made for a few nuclei, as shown in Figures 1-5. The 
photonuclear cross sections in this section are calculated using experimental absorption cross 
sections, from reference [19], at various photon energies. The experimental cross sections are 
multiplied by the WE branching ratio, given in equation (14), to produce the photonuclear 
cross sections shown in red in Figures 1-5. Similarly, the energy independent branching 
ratios are used, equations (20) and (21), to find the photonuclear cross section results shown 
in green. 

By using experimental absorption cross sections, the WE and energy independent branch- 
ing ratio can be compared directly, because it eliminates the possibility of induced uncer- 
tainty from a parameterized absorption cross section. The absorption cross section used in 
this paper, equation (2), for EMD nucleus-nucleus cross section comparisons is a Lorentzian 
parameterization. This absorption cross section will give the photonuclear cross section its 
dominant feature, a hump, which corresponds to the excitation of the giant dipole resonance. 
This parameterization will be poor for deformed nuclei, where the cross section is seen as 
the sum of two Lorentzians, and for light nuclei, which often exhibit complicated shapes 
[ 20 ]. 
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It can be seen in Figures 1-5 that the theoretical photonuclear cross sections, calculated 
with WE and energy independent branching ratios, adequately agree with experiment. For 
most of the reactions, the theoretical photonuclear cross section, calculated with a WE 
branching ratio, provides better agreement to experiment. This can be observed for the 
reactions 28 Si( 7 ,n), 28 Si( 7 ,p), 88 Sr( 7 ,n), and 91 Zr(y,n), as shown in Figures 2-5, respectively. 
It should be noted that for Figures 4 and 5, the WE branching ratio provides a much better 
fit to the experimental values around the peak of the photonuclear cross section, therefore 
capturing the essential feature of the GDR. 

For 16 0(7,n), shown in Figure 1, it is unclear which branching ratio (WE or energy 
independent) provides the better fit to experimental results, without performing a statistical 
analysis. The WE theory will have problems evaluating reactions involving oxygen, since a 
significant contribution of the photonuclear reaction comes from direct reactions [19]. The 
WE theory is not intended to be used for reactions that proceed directly and do not pass 
through a well defined compound state. 

Problems matching experimental data can be seen for reactions involving 28 Si. The 
theoretical photoneutron and photoproton cross sections calculated using WE branching 
ratios seem to provide the better fit. Although, as seen in Figures 2 and 3, the WE branching 
ratio does not represent the data well. Reactions involving 28 Si are unique in the fact that 
compound nucleus decay is found for the neutron channel, but different behavior is observed 
for proton decay. Predominately direct Giant Dipole Resonance (GDR) decay occurs for 
the proton channels of 28 Si [21]. Most of the partial (ground state, first excited state, etc.) 
proton cross sections demonstrate different intermediate structures [21]. Consequently, the 
branching ratios calculated with WE theory will not be representative of the proton decay 
channel. The erroneously calculated proton branching ratio will also have an effect on the 
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other channels (neutron, deuteron, triton, helion, and alpha), because the denominator of 
the WE branching ratio formula contains the total width, which is given by the sum of all 
the decay channels. If one of the channels is calculated incorrectly, it will ultimately affect 
all of the channels. 

From the comparisons presented, it is recommended that NUCFRG2 should use WE 
branching ratios rather than the energy independent method currently employed. For com- 
pound nucleus reactions, as in equation (13), the WE branching ratio provides the best 
fit, which is shown for 88 Sr( 7 ,n) and 91 Zr( 7 ,n). For reactions that have single or multiple 
channels that proceed through a direct reaction, as in equation (19), the WE branching 
ratio provides the same or better results when compared to the energy independent branch- 
ing ratio. This is surprising, since WE theory was not designed for reactions that proceed 
through direct channels. Due to the WE branching ratio giving similar, or better, results 
than the energy independent branching ratio, for direct and compound nucleus reactions, 
it is recommended that WE branching ratios be incorporated into NUCFRG2. In addition, 
the WE branching ratios are more advantageous, because they allow for not only neutron 
and proton emission, but also for deuteron, triton, helion, and alpha particle emission. 

B. Nucleus-Nucleus Cross Sections using WE and energy independent branching 
ratios 

In the previous section, a comparison of WE and energy independent branching ratios 
was presented for photonuclear cross sections. Conclusions made from the comparison and 
evaluation of photonuclear cross section results should be used for making future NUCFRG2 
recommendations. This comparison eliminated the possibility of induced uncertainty from 
parameterized terms, as mentioned previously. In this section, EMD cross sections for 
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single neutron, proton, and alpha particle removal are calculated using WE and energy 
independent branching ratios and then compared to experiment, as shown in Table III. 
Graphical depictions of this data are shown in Figures 6-20. 

To calculate nucleus-nucleus cross sections using an energy independent branching ratio, 
equation (24) is employed, where the proton and neutron branching ratios are given by 
equations (20) and (21), respectively. It is important to reiterate that for the case of the 
energy independent branching ratio, the formulation was based on the assumption that only 
neutron and proton emission occurred. Therefore, single alpha particle removal EMD cross 
sections, calculated using energy independent branching ratios, could not be determined. In 
Table III, for single particle alpha particle removal, the results calculated with an energy 
independent branching ratio will be marked not applicable (N/A). The current capability 
to calculate only neutron and proton emission in NUCFRG2, due to energy independent 
branching ratios, makes WE theory very desirable. 

For a nucleus-nucleus cross section that utilizes the WE branching ratio, equation (23) 
is used. The WE branching ratio is given by equation (14), with the total width given by 
equation (15) and the partial width given by equations (16) and (17) for a neutron and 
charged particle, respectively. Here, EMD cross sections for single alpha particle removal 
that employ WE branching ratios can be calculated, and are given in Table III. 

Before comparing EMD nucleus- nucleus cross section results to experiment, using WE 
and energy independent branching ratios, it is necessary to consider a few caveats in this 
analysis. One is that EMD nucleus-nucleus cross sections are calculated with a parameterized 
absorption cross section. This parameterization is not intended for deformed nuclei, such 
as 28 Si [22] or for light nuclei, which exhibit complicated absorption cross section shapes. 
Another caveat is that there is very little experimental data concerning EMD nucleus-nucleus 
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cross sections. Most experimental efforts have focused on special cases, in which the decay 
channel displays interesting behavior. Consequently, the experimental EMD nucleus-nucleus 
cross section data are not a good representation of a typical reaction. In addition, many 
of these experimental results have large errors, which demand better statistics for more 
accurate comparisons. It is important to consider all of these caveats when comparing EMD 
nucleus-nucleus cross section results, using WE and energy independent branching ratios, 
to experiment. For all the reasons listed above, the following results should not be used for 
making recommendations for NUCFGR2, but can offer insight into EMD nucleus- nucleus 
cross sections results that would be predicted by WE and energy independent branching 
ratios. 

It can be seen that the EMD cross sections with decay channels n C + In and n B + 
lp, for a 12 C projectile, calculated using WE and energy independent branching ratios, 
are within the experimental error estimates for mostly all target and energy combinations, 
as shown in Figures 6-9. This also occurs with 15 0 + In and 15 N + lp decay channels, 
for 16 0 projectiles, as shown in Figures 10 and 11. It is surprising that the EMD cross 
sections calculated with WE branching ratios agree so well with the experimental data 
because 12 C and 16 0 reactions proceed mainly through a direct reaction. The WE theory 
is specifically designed for compound nucleus decay and not intended for direct decay, as 
mentioned previously. Note that for the 15 N + lp decay channel for a 16 0 projectile has a 
negative experimental cross section, which is due to the large error bars. 

It is obvious that discrepancies exist between the calculated EMD cross sections and 
experimental data for 18 0 projectiles, as shown in Figures 12 and 13. The EMD cross sections 
calculated with WE and independent branching ratios will show disagreement because the 
parameterization of the absorption cross section, equation (2), was designed for only stable 
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nuclei. Before comparisons can be made for 18 0 projectiles, a parameterization for the 
absorption cross section for unstable nuclei will need to be developed. 

Unexpected behavior can be observed for the 28 Si projectile. The EMD cross sections 
are calculated with a parameterized absorption cross section. As mentioned previously, this 
parameterization is not intended for a deformed nuclei, such as 28 Si [22]. It is surprising 
that the results predicted with a WE or energy independent branching ratio, give reasonable 
results for single neutron, proton, and alpha particle decay, as shown in Figures 14-16. In 
addition, WE results will be affected by the predominately direct GDR decay in the proton 
channel [21]. Therefore, the proton channel will be calculated incorrectly; this will ultimately 
effect the rest of the channels, as mentioned in the previous section. 

Currently, NUCFRG2 does not have the capability of calculating alpha particle re- 
moval EMD cross sections. The energy independent branching ratio, currently employed 
in NUCFRG2, can calculate neutron and proton emission. In addition to neutron and 
proton removal, the WE branching ratio offers not only the calculation of alpha particle 
removal, but also deuteron, hclion, and triton removal. The EMD cross sections for alpha 
particle removal of 28 Si on various targets, displayed in Figure 16, show results that are 
predicted by WE theory. It would be advantageous for NUCFRG2 to calculate additional 
channels, such as the alpha channel. Therefore, it is important to take a closer look at 
the predicted EMD cross sections for alpha particle removal. In Figure 16, it is observed 
that the theoretical cross sections, calculated with a WE branching ratio, compare well to 
experiment for 2, A1, 64 Cu, and 120 Sn targets. Although, WE theory gives a much larger 
cross section than experimentally predicted for a 208 Pb target. It is proposed here that this 
experimental EMD cross section for single alpha particle removal of 28 Si is too small. The 
scattering cross section should be proportional to the square of the atomic number of the 
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target. Using the principle of detailed balance, a simple estimate of the EMD cross section 
for single alpha particle removal of 28 Si on a 208 Pb target can be calculated. From the rela- 
tion <7208 Pb oc ^ z ^ospb j w i iere <j 208 pf, and Z 2 ospb are the cross section and atomic number 

for a 208 Pb target, respectively. For <7* and Z*, the cross section and atomic number values 
of 2 'A1, 64 Cu, and 120 Sn targets are used (the cross section data are taken from Table 111). 
This simple and quick estimate gives a cross section value of approximately 836 mb, 136 
mb, and 156 mb when using 2 'A1, 64 Cu, and 120 Sn, respectively. The range of the estimated 
EMD alpha particle removal cross section of 28 Si on a 208 Pb target is 136-836 mb, which 
is higher than the experimental value of 72 ± 32 mb. It would be premature to invalidate 
the EMD alpha particle removal cross section of 28 Si on a 208 Pb target, calculated with WE 
theory (264.33 mb), because its value lies well within the estimated range of 136-836 mb. 
Until better statistics or more experimental data becomes available, accurate comparisons 
cannot be made. Consequently, there is no reason to discredit using WE theory to calculate 
alpha particle removal EMD cross sections. 

Problems calculating EMD cross sections with WE branching ratios can be observed 
for the 88 Y + In decay channel, Figure 17. This occurs for every projectile and energy 
combination. The energy independent branching ratio offers reasonable estimates, while the 
WE branching ratio experiences problems. Complications arise in 89 Y because 75% of the 
photoproton cross section between 14 and 24 MeV decays directly [23-25] . The partial direct 
decay behavior in the proton channel will not be accounted for by the WE theory, which 
calculates decay rates for compound nucleus reactions. Just as in the case of 28 Si projectiles, 
the erroneously calculated proton channel will affect all of the other decay channels. 

For the decay channels, 58 Co + In and 19 'Au + In, the WE and energy independent 
branching ratios give similar results for mostly all of the projectile and energy combinations, 
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Figures 18 and 19-20, respectively. It is important to note that out of all the EMD cross 
sections in this work, these two decay channels are the only channels that follow compound 
nucleus decay and do not have one or more channels that follow direct decay. This means that 
these are the only decay channels in this analysis that the WE theory is intended to apply. 
Both the WE and energy independent branching ratio calculate EMD cross sections, for the 
58 Co + In and 197 Au + In decay channels, that agree reasonably well with experiment. 

VII. CONCLUSION 

The nuclear fragmentation model NUCFRG2 currently calculates EMD nucleus- nucleus 
cross sections using energy independent branching ratios. Results indicate that the WE 
theory provides photonuclear cross sections that are comparable, but in most cases, better 
than the results obtained with energy independent branching ratios. Consequently, the WE 
theory should be considered as a good alternative to the energy independent branching ratio, 
which is currently employed in NUCFRG2. It was argued that the NUCFRG2 model should 
incorporate WE theory for the calculation of EMD cross sections and should consider the 
sensitivity of the results over computational time. The energy independent branching ratio 
allows for faster computation times of EMD nucleus-nucleus cross sections, due to its energy 
independence and simple formalism. The ability of the WE branching ratio to calculate 
not only neutron and proton emission, but also deuteron, triton, helion, and alpha particle 
emission make WE theory more advantageous. 

Future work will include the investigation of a branching ratio formalism that can be 
applied to direct channels for future implementation into NUCFRG2. 
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TABLE I: Values for the width of the electric dipole giant 
dipole resonance. 


Nucleus 

T (MeV) 

Reference 

12 C 

8.0 

[5] 

16 0 

10.0 

[5] 

18 0 

12.0 

[5] 

28 Si 

10.0 

[5] 

59 Co 

9.0 

[26] 

89y 

3.95 

[27] 

197 Au 

3.5 

[5] 
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TABLE II: Values for the 10% Charge Density Radius found 
from various models. 


Nucleus 

10% Chage Radius (frn) 

Reference 

12 C 

3.33 

[5] 

16 0 

3.77 

[5] 

18 0 

3.88 

[5] 

20 Ne 

4.06 

[5] 

27 A1 

4.21 

[5] 

28 Si 

4.18 

[5] 

32 S 

4.53 

[5] 

48 Ti 

5.00 

[5] 

56 Fe 

5.28 

[5] 

59 Co 

5.33 

[28] 

64 Cu 

5.45 

[5] 

89y 

6.02 

[28] 

108 Ag 

6.32 

[5] 

120 Sn 

6.58 

[28] 

139 La 

6.89 

[29] 

197 Au 

7.56 

[5] 

208pp 

7.83 

[5] 

2° 9 Bi 

7.78 

[28] 

238p- 

8.13 

[5] 
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TABLE III: Electromagnetic Cross Sections for Single Neu- 
tron, Proton, and Alpha Particle Removal, are the 

experimental EMD cross sections from references [30-36]. 

are the theoretical EMD cross sections using the en- 
ergy independent branching ratio, equation (24). ct™ d are 
the theoretical EMD cross sections using a WE branching 
ratio, equation (23). Note that all data from this table are 
graphically displayed in Figures 6-20. 


Projectile Target Tj a b (GeV/N) Decay Channel <r™ t D (mb) cr™ D (mb) cr||^ D (mb) 


12 c 

208 pb 

2.1 

n C + In 

51 ± 18 

42.9964 

35.2666 

12 c 

208 pb 

2.1 

n B + lp 

50 ± 25 

70.2699 

38.6367 

12 c 

208 pp 

1.05 

n C + In 

39 ± 24 

25.2166 

19.8329 

12 c 

208 pp 

1.05 

n B + lp 

50 ± 25 

41.8808 

22.1407 

16 o 

208 pp 

2.1 

15 0 + In 

50 ± 24 

62.051 

50.2401 

16 0 

208 pb 

2.1 

15 N + lp 

96 ± 26 

102.91 

83.5714 

12 c 

108 Ag 

2.1 

n C + In 

21 ± 10 

16.7542 

13.8676 

12 c 

108 Ag 

2.1 

n B + lp 

18 ± 13 

27.2921 

15.1292 

12 c 

108 A g 

1.05 

n C + In 

21 ± 10 

10.5761 

8.44344 

12 c 

108 Ag 

1.05 

n B + lp 

25 ± 19 

17.4608 

9.36414 

16 0 

108 Ag 

2.1 

15 0 + In 

26 ± 13 

23.9174 

19.532 

16 0 

108 Ag 

2.1 

15 N + lp 

30 ± 16 

39.4881 

32.1429 


12 C 

64 Cu 

2.1 

n C + In 

10 ± 7 

7.06394 

5.87764 

12 C 

64 Cu 

2.1 

n B + lp 

4 ± 8 

11.4855 

6.39671 

12 C 

64 Cu 

1.05 

n C + In 

9 ± 8 

4.65276 

3.74663 

12 C 

64 Cu 

1.05 

n B + lp 

5 ± 8 

7.65609 

4.13926 

16 0 

64 Cu 

2.1 

15 0 + In 

9 ± 8 

10.0187 

8.2224 

16 0 

64 Cu 

2.1 

15 N + lp 

15 ± 8 

16.499 

13.4477 

12 c 

27 A1 

2.1 

14 C + In 

0 ± 5 

1.6516 

1.3846 

12 c 

27 ai 

2.1 

n B + lp 

0 ± 5 

2.67839 

1.50157 

12 c 

27 A1 

1.05 

14 C + In 

1 ± 6 

1.15605 

0.942364 

12 c 

27 A1 

1.05 

n B + lp 

1 ± 7 

1.89353 

1.03538 

16 o 

27 A1 

2.1 

15 0 + In 

0 ± 5 

2.31977 

1.91739 

16 0 

27 A1 

2.1 

15 N + lp 

-1 ± 9 

3.80662 

3.10831 


Continued on Next Page. . . 
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TABLE III - Continued 


Projectile Target Tj a b (GeV/N) Decay Channel 

(mb) 

u™ (mb) 

(mb) 

12 C 

12 c 

2.1 

n C + In 

-2 ± 5 

0.393841 

0.331945 

12 C 

12 c 

2.1 

n B + lp 

-1 ± 4 

0.637513 

0.359076 

12 C 

12 c 

1.05 

n C + In 

-2 ± 5 

0.287878 

0.236711 

12 C 

12 c 

1.05 

n B + lp 

-2 ± 5 

0.470029 

0.259044 

16 0 

12 c 

2.1 

15 0 + In 

-1 ± 4 

0.549095 

0.456151 

16 0 

12 c 

2.1 

15 N + lp 

-1 ± 4 

0.898782 

0.734837 


18 0 

238 pj 

1.7 

17 0 + In 

140.8 ± 4.1 

87.5476 

211.602 

18 0 

238-y 

1.7 

17 N + lp 

25.1 ± 1.6 

101.782 

0.419803 

18 0 

208pp 

1.7 

17 0 + In 

136 ± 2.9 

71.7717 

173.512 

18 0 

208pp 

1.7 

17 N + lp 

20.2 ± 1.8 

83.6706 

0.34833 

18 0 

48 Ti 

1.7 

17 0 + In 

8.7 ± 2.7 

7.05503 

17.0903 

18 0 

48 Ti 

1.7 

17 N + lp 

-0.5 ± 1.0 

8.43365 

0.038370S 


28 Si 

208pp 

14.6 

27 Si + In 

241.0 ± 4.2 

409.25 

153.513 





263.5 ± 15 



28 Si 

208 Pb 

14.6 

27 A1 + lp 

676.4 ± 7.6 
672.2 ± 25 

516.844 

646.874 

28 Si 

208pp 

14.6 

24 Mg + la 

72 ± 32 

N/A 

264.33 


28 Si 

120 Sn 

14.6 

27 Si + In 

100.0 ± 2.0 
101.9 ± 8.2 

160.077 

60.1917 

28 Si 

120 Sn 

14.6 

27 A1 + lp 

274.0 ± 4.4 

201.801 

252.733 





290.6 ± 13 



28 Si 

120 Sn 

14.6 

24 Mg + la 

58 ± 30 

N/A 

102.647 


28 Si 

64 Cu 

14.6 

27 Si + In 

40.4 ± 1.7 

56.5488 

21.3114 

28 Si 

64 Cu 

14.6 

27 A1 + lp 

111.0 ± 3.3 

71.17 

89.1858 

28 Si 

64 Cu 

14.6 

24 Mg + la 

17 ± 19 

N/A 

36.0171 


28 Si 

27 A1 

14.6 

27 Si + In 

13.11 ± 0.59 

12.0398 

4.54909 





10.8 ± 2.8 



28 Si 

27 A1 

14.6 

27 A1 + lp 

31.6 ± 1.2 
34.2 ± 4.7 

15.1243 

18.9656 

28 Si 

27 A1 

14.6 

24 Mg + la 

21 ± 16 

N/A 

7.6096 


Continued on Next Page. . . 
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TABLE III - Continued 


Projectile Target Tj a b (GeV/N) Decay Channel a™ t D (mb) cr™ D (mb) cr||^ D (mb) 


12 C 

197 Au 

2.1 

196 Au + In 

75 ± 4 

36.9002 

34.2595 

12 C 

89y 

2.1 

88 Y + In 

9 ± 12 

11.9561 

0.152798 

12C 

59 Co 

2.1 

58 Co + In 

6 ± 9 

5.37801 

3.21086 

16 0 

197 Au 

60 

196 Au + In 

280 ± 30 

214.465 

201.668 

16 0 

197 Au 

200 

196 Au + In 

440 ± 40 

278.427 

262.135 

20 Ne 

197 Au 

1.7 

196 Au + In 

151 ± 13 

97.2338 

90.2015 

20 Ne 

197 Au 

2.1 

196 Au + In 

153 ± 18 

88.6733 

82.084 

20 Ne 

89y 

2.1 

88 Y + In 

43 ± 12 

31.1781 

0.392326 

20 Ne 

59 Co 

2.1 

58 Co + In 

32 ± 11 

13.9518 

8.3017 

40 Ar 

197 Au 

1.8 

196 Au + In 

348 ± 34 

279.913 

259.048 

40 Ar 

89y 

1.8 

88 Y + In 

132 ± 17 

88.5417 

1.06924 

56 Fe 

197 Au 

1.7 

196 Au + In 

601 ± 54 

545.735 

504.356 

56 Fe 

89y 

1.7 

88 Y + In 

217 ± 20 

170.922 

2.01472 

56 Fe 

59 Co 

1.7 

58 Co + In 

88 ± 14 

75.779 

44.447 

139 La 

197 Au 

0.15 

196 Au + In 

447 ± 28 

604.729 

510.793 

139 La 

197 Au 

1.26 

196 Au + In 

1970 ± 130 

1994.08 

1830.42 

139 La 

59 Co 

1.26 

58 Co + In 

280 ± 40 

263.184 

150.123 

32 S 

197 Au 

200 

196 Au + In 

1120 ± 160 

1099.77 

1035.37 

197 Au 

197 Au 

1.0 

196 Au + In 

3077 ± 200 

3207.16 

2927.37 

209 Bi 

197 Au 

1.0 

196 Au + In 

3244 ± 205 

3469.23 

3165.08 

238-g 

197 Au 

0.96 

196 Au + In 

3160 ± 230 

4037.1 

3677.04 




_i_ i ^ : i i i i i i i i i i 

15 20 25 30 

E y (MeV) 


FIG. 1: Comparison of theoretical and experimental photoneutron cross sections for 16 0. The 
theoretical photoneutron cross sections are calculated using WE (red) and energy independent 
branching ratios (green). Experimental data are from Figure 14(c) of reference [27]. 
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28 Si(y,n) 



FIG. 2: Same as Figure 1, except for 28 Si. Experimental data are from Figure 15(b) of reference 
[27]. 


30 


28 si(r,p) 



FIG. 3: Same as Figure 1, except for the photoproton cross section for 28 Si. Experimental data 
are from reference [19] (p. 110). 
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88 Sr(y,n) 



32 


91 Zr(y,n) 



FIG. 5: Same as Figure 1, except for 91 Zr. Experimental data are from reference [19] (p. 171). 
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Projectile: 12 C, 7j ab = 2.1 (GeY/N), Decay Channel: n C + In 
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-10 

12 C 27 A1 64 Cu 108 Ag 208 Pb 
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FIG. 6: Comparison between theory and experiment for the single neutron removal EMD cross 
section of 12 C for various targets at 2.1 GeV/N. Experimental data [displayed as black circles] are 
from reference [32]. Theoretical cross sections are calculated using WE [displayed as red squares, 
using equation (23)] and energy independent [displayed as green triangles, using equation (24)] 
branching ratios. Note that this is a graphical representation of the data in Table III. 
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Projectile: 12 C, 7j ab = 2.1 (GeV/N), Decay Channel: U B + lp 
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FIG. 7: Same as Figure 6, except for single proton removal. 
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Projectile: 12 C, 7j ab = 1.05 (GeV/N), Decay Channel: n C + In 
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FIG. 8: Same as Figure 6, except at 1.05 GeV/N. 
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Projectile: 12 C, 7j ab = 1.05 (GeV/N), Decay Channel: n B + lp 



12 C 27 A1 64 Cu 108 Ag 208 Pb 


Target 


FIG. 9: Same as Figure 6, except for single proton removal at 1.05 GeV/N. 
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Projectile: 16 0, 7j ab = 2.1 (GeV/N), Decay Channel: 15 0 + In 
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• cr(E y ) Experiment 
o-(E y ) = g n JdE y N(E y ) cr abs (E y ) 

□ o-(Ey) = /dE y N (Ey) g n (E ) o- ahs (E ) 
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FIG. 10: Same as Figure 6, except for 16 0 at 2.1 GeV/N. 
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Projectile: 16 0, 7j ab = 2.1 (GeV/N), Decay Channel: 15 N + lp 
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• cr(E y ) Experiment 
o-(E y ) = g p f dE y N(E y ) cr abs (E y ) 

□ <r(Ey ) = fdEy N(Ey ) g p (Ey ) o- abs (Ey ) 



FIG. 11: Same as Figure 6, except for the single proton removal of 16 0 at 2.1 GeV/N. 
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Projectile: ls O, 7j ab = 1.7 (GeV/N), Decay Channel: n O + In 
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FIG. 12: Same as Figure 6, except for 18 0 at 1.7 GeV/N. Experimental data [displayed as black 
circles] are from reference [34]. 
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Projectile: 18 0, 7j ab = 1.7 (GeV/N), Decay Channel: 17 N + lp 
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FIG. 13: Same as Figure 6, except for the single proton removal of 18 0 at 1.7 GeV /N. Experimental 
data [displayed as black circles] are from reference [34] . 
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Projectile: 28 Si, 7j ab = 14.6 (GeV/N), Decay Channel: 27 Si + In 
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FIG. 14: Same as Figure 6, except for 28 Si at 14.6 GeV/N. Experimental data [displayed as black 
circles] are from references [30, 31, 37]. 
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Projectile: 28 Si, 7j ab = 14.6 (GeV/N), Decay Channel: 27 A1 + lp 
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FIG. 15: Same as Figure 6, except for single proton decay of 28 Si at 14.6 GeV/N. Experimental 
data [displayed as black circles] are from references [30, 31, 37]. 
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Projectile: 28 Si, 7j ab = 14.6 (GeV/N), Decay Channel: 24 Mg + la 
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FIG. 16: Same as Figure 6, except for single alpha particle decay of 28 Si at 14.6 GeV /N. Experi- 
mental data [displayed as black circles] are from reference [30]. Note that there are no theoretical 
cross sections calculated with an energy independent branching ratio since an energy independent 
branching ratio has not been formulated for alpha particle removal. 
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Target: 89 Y, Decay Channel: 88 Y + In, Projectiles: 12 C, 2C) Ne, 40 Ar, and 56 Fe 



Projectile at 7j ab (GeV/N) 

FIG. 17: Comparison between theory and experiment for the single neutron removal EMD cross 
section of a 89 Y target for various projectile and energy combinations. Note that this is a graphical 
representation of the data in Table III. Experimental data [displayed as black circles] are from 
reference [35]. Theoretical cross sections are calculated using WE [displayed as red squares, using 
equation (23)] and energy independent [displayed as green triangles, using equation (24)] branching 
ratios. 
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Target: 59 Co, Decay Channel: 58 Co + In, Projectiles: 12 C, 20 Ne, 56 Fe, and 139 La 
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FIG. 18: Same as Figure 17, except for a 59 Co target for various projectile and energy combinations. 
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Target: l97 Au, Decay Channel: 196 Au + In, Projectiles: 12 C, 16 0, and 20 Ne 
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FIG. 19: Same as Figure 17, except for a 197 Au target for various projectile and energy combina- 
tions. Experimental data [displayed as black circles] are from references [36] and [35]. 
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Target: 197 Au, Decay Channel: 196 Au + In, Projectiles: 32 S, 40 Ar, 56 Fe, and 139 La 
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FIG. 20: Same as Figure 17, except for a 19 'Au target for various projectile and energy combina- 
tions. Experimental data [displayed as black circles] are from references [36] and [35]. 
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